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■ Suppose 



1. Introduction, basic idea and results. 

A(u) = f, 



(1.1) 
(1.2) 



O 

(N : a(o) = o, 

where A(u) is a nonlinear operator on a Banach space X , which has the G- derivative: 
Q : A'{u)w = lim + yw 

o 



Here D(A) is the domain of definition of the operator A, which is assumed to be a set, such that if u 
and w belong to it, then for all sufficiently small t £ M., the element u + tw belongs to it. The operator 
A'(u) is a linear, maybe unbounded, densely defined operator on X. Even when A'(u) is unbounded, it 
may be that the following estimate holds: 



o 
o 

43 ; \\A'(u) - -> as \\u - v\\ -> 0. 

This is usually the case when the unbounded part of A'(u) is an operator independent of u, see Example 
2 in section 2 below. 

In this paper we develop an idea from [1] to write (1.1) as an equation: 

U 1 L(u)u = /, A{u) = L(u)u, (1.4) 



where L{u) is a linear operator on X which depends nonlinearly on u(x) as on a parameter. We assume 
that L(u) has bounded inverse, that is L^ 1 (m) is a bounded linear operator defined on all of X. 

We study an approach described in [1] for solving equation (1.1) written in the form (1.4). This 
approach differs from the traditional ones (such as Newton's method and its modifications, and other 
methods, based on local linearizations of nonlinear operator A{u) (sec [2] for a detailed study of such 
methods). The difference is in the global nature of equation (1.4) with linear operator L{u) which depends 
nonlinearly on u. 

In [1] the goal was to construct an analog of Green's function for nonlinear systems and to suggest 
numerical schemes, which would be based on the schemes for solving linear equations. The hope is 
that such schemes may be more efficient and more stable computationally. In [1] an iterative scheme 
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(1 . 1 1)- (1 . 12) with uo = f was suggested. Analysis of the convergence of this iterative scheme was not 
given in [1] but numerical examples in [1] are encouraging. Here wc give such an analysis. 

In Theorem 1.1 we give sufficient conditions for convergence of an iterative process (1.11)-(1.12) (see 
below) for solving equation (1.1). These conditions are formulated in terms of the operator L(u) defined 
in (1.6) below. 

In Theorem 1.2 sufficient conditions are given in terms of the G-derivative of A(u) for the operator 
L(u) to satisfy the conditions of Theorem 1.1. 
In section 2 examples are briefly considered. 

For linear equation (1.4) one can use various methods for constructing L _1 (m). One writes (1.4) as 

u = L- 1 (u)f, (1.5) 

assuming that L(u) is boundedly invertible ( sufficient conditions for this are given in Theorem 1.2) and 
uses an iterative scheme to solve equation (1.5). 
Let us first prove that 

L{u) = ( A'(tu)dt, (1.6) 
Jo 

and that 

A{u) = I A\tu)dtu = L{u)u. (1.7) 
Jo 

The following lemma is an immediate consequence of the definition of the G-derivative and its proof is 
included for convenience of the reader. 

Lemma 1.1. Formula (1.7) holds. 

Proof of Lemma 1.1. Note that 

/ A(tu) dt = A(u) - A(0) = A(u), (1.8) 
Jo dt 

where the assumption (1.2) was used. 
On the other hand, 

d . . , ,. A((t + s)u) -A(tu) 

~ A(tu) = lim — ^ — = A'{tu)u, (1.9) 

dt s^o s 

by the definition of the G-derivative. It follows from (1.8), (1.9), and the linearity of A'(u) as an operator 
on X, that 

ni 

A'{tu)dtu = A(u), (1.10) 
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f 

Jo 



as claimed. 

Lemma 1.1 is proved. □ 

Now the linear operator L(u) on X (depending on u nonlincarly) is well defined by formula (1.6), 
and equation (1.4) is well understood as an equation for u with linear operator L(u), which depends 
nonlinearly on u. 

If L(u) is boundedly invertible, that is, _L _1 (u) is a bounded linear operator defined on all of X, then 
(1.4) is equivalent to (1.5). Below we use the boundedness of L _1 (m) not for all u G X but only in a ball 
Br ■= \\u\\ < R, where R > is some number which will be specified in the proof of Theorem 1.1. 

Let us write an iterative process for solving equation (1.5): 

U n +1 = L _1 (u„)/, (1.11) 



MO = MO, 
2 



(1.12) 



where uq is an initial approximation which is arbitrary at this moment. In practice it should be chosen 
taking into account all the a priori information available about the solution. This may improve the rate 
of convergence of the scheme (1.11)-(1.12). 

Let us now analyze the above iterative process. 

We have 

||u n+ i - u n \\ = HL" 1 ^) - L _1 (m„_i)|| 11/11 < q\\u n - u„_i|| ||/||. (1.13) 



Here we have assumed that: 



|L _1 (u) — L _1 (w)|| < q\\u- v\\, 



(1.14) 



for u, v G Br, where Br is a ball < R, R > is a certain number, and g > is a number independent 
of u, v G _B_r. 

Assume now that 

Then 



<Q < 1. 



(1.15) 



||«i-«o|| < ||£ _1 («o)/-«o||, 
||«2 - "ill < - £ _1 («o)|| 11/11 < «||/|| ||«i - MoH < QHui - Uo| 

U n \\ < WL-^Un) - L _1 (u„_i)|| ll/H < Q\\u n - U„_l|| 

< <3 2 ||u„-i - u„_2|| < • • • < Q"||«i - "oil- 



\Un+l 



Thus 



l«n+l|| = 



n+1 



+ - u 3-l] 



^ II II , II II i-Q"^ ^ || || . <m 
< ||u || + IK - "oil — : 7T~ ^ IKII + 



(1.16) 



1-Q " 11 u " 1-Q 

It is sufficient to assume that conditions (1.14) -(1.15) hold in a ball Br such that Br d S, where 

S:={u: \\u\\ < \\u \\ + (1 - Q)" 1 ^! - u ||}. (1.17) 

If Br D S, then all the members of the sequence u n belong to S, and the sequence u n converges to the 
element 



u := uq + ^~^(uj — Uj-i), u G S. 

3 = 1 



(1.18) 



Passing to the limit in (1.11) and using condition (1.14), which implies the continuity (in the operator 
norm) of L _1 (u) with respect to u, one gets 



u = L- 1 {u)f, 



(1.19) 



or 

A(u) = L{u)u = f. (1.20) 

Let us formulate the result. 

Theorem 1.1. Assume that A(u) has the G-derivative A'(u) which is a linear, possibly unbounded, 
densely defined operator on X, such that for the operator L(u), defined in equation (1-6), the conditions 
(1-14) o,nd (1.15) hold for allu and v in S, where S is defined in (1.17). Then equation (1.1) is uniquely 
solvable in S and its solution can be obtained by the iterative process (1.11)-(1.12). 

Remark 1.1. One could choose uq = f in (1.12). However, in practice there may be some information 
available which allows one to choose a better initial approximation uq(x) to the unknown solution u. 
Note that our conditions (1.14)-(1.15) are the usual conditions for the applicability of the contraction 
mapping principle to the equation u = i _1 (u)/. 
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One could formulate a condition, which implies (1.15), in terms of the bounds on the derivative of 
L _1 (u) with respect to u. Namely L~ 1 (u) u = — L _1 (m)L(u) u _L _1 (m), where the subscript u denotes the 
G-derivative. Therefore (1.14) holds with q = sup u ||L _1 (ti)L(u) u L _1 (it)||. We do not go into further 
detail since conditions (1.14), (1.15) are simple and can be verified in some examples (see for instance 
example 2 in section 2). 

Note that an application of the contraction mapping principle directly to the equation (1.1) is not 
possible in many practical problems in which the operator A(u) is unbounded. Therefore equation (1.5) 
is much more suitable for an application of the contraction mapping principle in these problems. 

Remark 1.2. The inequality, which implies that S C Br, is of the form: 

IM+ i-«z||/|| (L21) 

Solving it for q, one gets the following inequality: 

.<.< 11/11- (,- M^pt ). a,,, 

For a given / (/ ^ since otherwise u = by the assumption (1.2) and the local injectivity of A(u) 
near u = 0, which follows from the assumed bounded invcrtibility of L(0) ), one can fix an arbitrary uo 
and then choose R, R > \\uo\\, sufficiently large, so that 

\\L-H«o)f-uo\\ <lt (L23) 



R — || wo 

and 



0<Q<1, Q:=l- llL ~y°]f- U4 . (1.24) 

R - \\uo\\ 



Then there exists q > which satisfies (1-22) and therefore (1.15) holds. 
In particular, if 

m = 0, and R > ||L _1 (0)/||, 



then equation (1.1) is solvable in the ball Br provided that (1.14) holds and 

o<9< n/ir 



Assumptions (1.14)-(1.15) imply the unique solvability of equation (1.1) in the ball Br. This follows 
from Theorem 1.1, but we give an independent simple proof of this claim. 
Indeed, suppose 

A{u)=L{u)u = f, A(v) = L(v)v = f, u,v e B R . (1.26) 

Then 

L(u)u = L(v)v, u = (L _1 («) - L-^v)) f + v, 
so, using (1.14) and (1.15), one gets 

||u - v|| < HL- 1 ^) - L-\v)\\ ll/H < Q \\u - v\\ < \\u - v\\. (1.27) 

It follows from (1.27) that u = v. This proves uniqueness of the solution to (1.1) in the ball Br. 
Existence of the solution follows from Theorem 1.1. 

Now we give sufficient conditions for the existence of the bounded inverse operator i _1 (u). 
We start with 
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Lemma 1.2. Assume that: 



i) [A'iO)]- 1 



<P, 



ii) sup 0<t<1 || A' (to) — A' (0)|| < s, Vu e Br := {u : \\u\\ < R}, where R > is some number, 
and 



Hi) ps < 1. 
Then 



sup 

0<t<l 



[A'itu)}- 1 



< 



1 — ps 



Vw e B R . 



(1.28) 



Proof. First note that in (1.28) existence of the bounded linear operator [A' (to)] 1 is also claimed. One 
has 

A' (to) = A'(0) + A' (to) - A'(0) = A'(0) \l + (A'(0)) _1 (yl'(to) - A'(0))j := A'(0)(I + B). 
Therefore 

[A l {tu)]- 1 = {I + B)-'[A'{Q)]-\ (1.29) 
where both inverse operators on the right-hand side of equation (1.29) exist, and by i) and ii), one has: 

\\B\\ <ps<l. 

Therefore: 



\\(I + B)-'\\< 



(1.30) 



I + (A'(0)) _1 / (A'(tu) - A'(0j) dt 
Jo 



1 — ps 

From (1.29), (1.30) and the assumption i), the desired conclusion (1.28) follows. Lemma 1.2 is proved. □ 

Now we want to prove that, under the assumptions of lemma 1.2, the operator L _1 (w) does exist and 
is bounded. 
Since 

L(u) = / A'(tu)dt, 
Jo 

one has 

L{u) = A'(0) 

Using i) and ii) of Lemma 1.2, and denoting [A'(0)] _1 := T, one gets 

L^iu) = [I + TM]- 1 T, 

where 

M:=f [A'(tu)-A'(u)]dt, \\M\\ < s, \\TM\\<ps<\. 
Jo 

Therefore 

IKJ + TM)- 1 !! < . (1.34) 

We have proved the following: 
Theorem 1.2. Assume i), ii) and Hi) of Lemma 1.2. Then L _1 (u) exists for all u G Br and 



(1.31) 

(1.32) 
(1.33) 



\\L-\u)\\ < 



1 — ps' 



(1.35) 



2. Examples. 



Example 1. 

Let X = C(D), the usual Banach space of continuous in the domain D functions, K(x,y,u) and 
K u (x,y,u) be continuous functions on D x D x R, and 

A(u) := u + B{u) :=«,+ / K (x, y, u(y)) dy, K(x, y, 0) = 0. 

Jd 

We have 

A'(u)w = w(x) + / — — (x,y,u(y))w(y)dy, 



D 



du 



and 



r r 1 gx f K (x v u(v)) 

L(u)w = w{x)+ J^j^ —{x,y,tu(y))dtw(y)dy = w(x) + y, ^ y " w{y)dy. 



This L{u) is a continuous linear operator on C(D) which depends on u(y) nonlinearly. 
Equation 

A(u) = f (2.1) 

is equivalent to 

L(u)u - /, (2.2) 

or, if L _1 (u) is a bounded linear operator, to 

u = L- 1 (u)f. (2.3) 



An iterative scheme 

u n +i = i _1 (u„)/, or L(u n )u n+1 = f, u = M , (2.4) 

can be studied and, under suitable assumptions on L(u), this scheme converges. 

Let us briefly outline a proof of the boundedncss of L _1 (u) in this example. Assume that K(x, y, u) = 
k(x,y)g(u), where k(x,y) is a continuous selfadjoint kernel, and g(u) is a continuous function such that 
> 0. Assume that u(x) > is a continuous function in the closure of D. Then L(u) is a Fredholm- 
type operator with index zero, and its null-space is trivial, as one can easily check using the above 
assumptions. By Fredholm's alternative, L _1 (m) is bounded. If the selfadjoint kernel k(x,y) has positive 
eigenvalues, and if g'(u) > m > 0, then equation (2.1) has no more than one solution. Indeed, if 
A(u) - A(v) = 0, then 

(u - v,g(u) - g(v)) + (k[g(u) - g(v)],g(u) - g(v)) = 

where the parentheses denote the inner product in L 2 (D), kw :— j D k(x, y)w(y)dy, and our assumptions 
imply that both terms in the above equations are nonnegative. Therefore they both vanish, and this 
implies that u = v as claimed. 

Example 2. 
Let 

A(u) := —Aw + g(x, u) = f(x), in D, u| g = 0, g(x,0) = 0, (2.5) 

where D C R n is a bounded domain with sufficiently smooth boundary S, f(x) is a given function, 
g(x, u) is continuous with respect to x e D and C 1 with respect to u e R. As A we take H = L 2 (D). 
One has 

A'(u)w = —Aw + g' u (x, u)w, 
6 



and 



L(u)w = —Aw + 



g(x,u(x))w 
u(x) 



(2.5) 



The operator L(u) here is unbounded and his unbounded part, —A, does not depend on u(x). Under 
suitable assumptions the operator L(u) is boundcdly invertible and one can estimate For 
example, suppose g(x, u) > and g ^"^ = a + q(x), where a = const > 0, and q(x) is a nonnegative 
function which depends on u(x). In this case 

L-\u) = (-A + a + qix))- 1 , HL- 1 ^)!! < 1. 

Here the operator —A + a + q(x) is defined in L 2 (D) by the Dirichlct boundary condition, or by its 
domain H 2 {D) n H X (D), where H m {D) are the Sobolev spaces. 
One can get the following estimate: 



\L-\u)-L-\v)\\ < \\L-\u)\\ \\L(u)-L(v)\\ \\L~\v) 



In (2.7) we have assumed that 



< 



c(q) = max 



g(u) g(y) 



c(q) I, 



< CO, 



and have used the following well-known identity: 



(2.7) 



(2.8) 



(2.9) 



with B\ = L(u), £?2 = L(v). Thus, the conditions (1.14)-(1.15), sufficient for the convergence of the 
iterative process (2.4), hold if c(g)||/||a~ 2 < 1. 

Example 3. 

Let D C R n be a bounded domain with a sufficiently smooth boundary S. Consider the problem 



du 



V • [a(u)Vu] = f(x,t) in D x [0,T], 
u \t=o = M o(a;), 



0. 



Write this problem as 

A(u) :— u(x, t) — I B(u) dr — u a (x) + / f(x,t)dr 
Jo Jo 

with u = u{x, t) e C 1 (h 2 (D), [0, T]j , and 

B(u) := V • [o(«)V«]. 



Assume that 



aeC 2 (R), 0<c<a<m, \a'\,\a"\<m, 



(2.10) 

(2.11) 
(2.12) 



(2.13) 

(2.14) 
(2.15) 



where c and m are some positive constants and the primes denote derivatives with respect to u. 
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One checks that 



where 
Denote 



A'(u)w = w — / B'(u)wdT, 
Jo 

B'(u)w = V • [a'(u)Vw] + a'{u){V 2 u)w + a"{u){Vufw + a'(u)Vu ■ Vw. 

r 

-f(u) := / a(v) dv. 
Jo 



Now one checks that 



[ B'(au)daw = V 
Jo 

' a{u) j(u) 



a(u) — a(0) 



Vw 



a(u) 7(m) 



u 



+ 



u 



Vu-Vw + 



a'(u) 2a(u) 2j{u) 



(Au) w 
(Vu) 2 w. 



Therefore, 



L(u)w = / A' (au) da w = w 



w 



o Jo 



a(u) — a(0) 



B'(au) da w dr 
a{u) 7(1*) 



u u z 



a{u) ~f(u) 
u u 2 



Vu-Vw + 



a'(u) 2a(u) 2~f{u 



(Au) w 

(Vu) 2 w S> dr 



(2.16) 
(2.17) 
(2.18) 



(2.19) 



(2.20) 



and this equation defines a linear unbounded operator L(u), which depends on u = u(x,t) nonlinearly 
as on a parameter. 

It seems at first sight that as u — > one has difficulties in the definition (2.20), since u is in the 
denominator. However, in fact there arc no difficulties: one checks using L'Hospital rule, for example, 
that the functions 

a{u) 7 (u) _ ua(u) - jju) _^ a'(0) ^ ^ 



u 



and 



are bounded as u — > 0. 



a'(u)u 2 - 2ua(u) + 2"/(u) 



o"(0) 



m^o 3 



(2.22) 
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